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ABSTRACT 

We present the first systematic study of (non-radio-sclcctcd) radio-loud 
narrow-line Seyfert 1 (NLSl) galaxies. Cross-correlation of the 'Catalogue of 
Quasars and Active Nuclei' with several radio and optical catalogues led to the 
identification of ~11 radio- loud NLSl candidates including 4 previously known 
ones. This study almost triples the number of known radio-loud NLSl galax- 
ies if all candidates are confirmed. Most of the radio-loud NLSls are compact, 
steep spectrum sources accreting close to, or above, the Eddington limit. The 
radio-loud NLSls of our sample are remarkable in that they occupy a previously 
rarely populated regime in NLSl multi- wavelength parameter space. While their 
[OIII]/H/? and Fell/H/? intensity ratios almost cover the whole range observed in 
NLSl galaxies, their radio properties extend the range of radio-loud objects to 
those with small widths of the broad Balmer lines. Their black hole masses are 
generally at the upper observed end among NLSl galaxies, but are still unusu- 
ally small in view of the radio loudness of the sources. Among the radio-detected 
NLSl galaxies, the radio index R distributes quite smoothly up to the critical 
value of -R ~ 10 and covers about 4 orders of magnitude in total. Statistics show 
that ~7% of the NLSl galaxies are formally radio-loud while only 2.5% exceed 
a radio index R > 100. Implications for NLSl models are discussed. Several 
mechanisms are considered as explanations for the radio loudness of the NLSl 
galaxies and for the lower frequency of radio-louds among NLSls than quasars. 
While properties of most sources (with 2-3 exceptions) generally do not favor 
relativistic beaming, the combination of accretion mode and spin may explain 
the observations. 

Subject headings: quasars: individual (TEXlllll+329, SDSSJ094857.3+002225, 
SDSSJ172206.03+565451.6, RXJ0134-4258, IRAS09426+1929, SBS1517+520, RXJ16290+4007, 
PKS0558-504, 2E1346+2637, IRAS20181-2244, RXJ23149+2243) - quasars: emis- 
sion lines - quasars: general - X-rays: galaxies - radio continuum: galaxies 
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1. Introduction 

The radio-loud radio-quiet bimodality of quasars is one of the long-standing unsolved 
problems in the research of Active Galactic Nuclei (AGN). While the presence of two classes 
of quasars, radio-louds and radio-quiets with a deficiency of sources at intermediate radio 
powers, was historically apparently well established (e.g., Kellerman et al. 1989, Miller et 
al. 1990, Visnovsky et al. 1992), other studies have questioned the mere existence of a radio 
dichotomy (e.g.. White et al. 2000, Cirasuolo et al. 2003), or have pointed out that larger 
samples are needed to put the bimodality on a firmer statistical basis (Hooper et al. 1995). 
On the other hand it was argued that carefully selected small, well-studied samples still 
reveal the bimodality (e.g., Sulentic et al. 2003), and that it is still present in a large sample 
of SDSS-FIRST galaxies (Iveciz et al. 2002). Approximately 15% of the quasar population 
is radio loud (Urry & Padovani 1995) for reasons that are still unknown. 

A second question that is presently under close scrutiny is whether or not there are 
relations between black hole (BH) mass and radio loudness, and especially, whether or not 
there is a limiting BH mass above which objects are preferentially radio loud, and whether 
or not radio-louds and radio-quiets show the same spread in their black hole masses (e.g., 
Laor 2000, Lacy et al. 2001, Oshlack et al. 2002, Woo and Urry 2002, Shields et al. 2003, 
McLure & Jarvis 2004, Metcalf & Magliocchetti 2006, Liu et al. 2006). 

While most attention so far concentrated on the radio-loud radio-quiet distinction of 
quasars, extending the studies to other subtypes of AGN may shed new light on both, the 
existence of a radio-loud radio-quiet bimodality and on the cause of radio loudness. We 
have started an investigation of the radio properties of Narrow-line Seyfert 1 galaxies. Apart 
from the topic just mentioned, the analysis of their radio properties is also of interest for 
addressing open questions in the study of NLSls, as detailed below. 

NLSl galaxies were identified by Osterbrock & Pogge (1985) as objects with small widths 
of the broad Balmer hnes. These go hand in hand with weak [OIII]5007/H/3toti emission and 
strong emission from Fell complexes (e.g., Boroson & Green 1992). Physical drivers of the 
NLSl phenomenon and correlations among emission-line and continuum properties are not 
yet well understood. There is growing evidence, though, that most NLSls are objects with 
high accretion rates, close to or even super-Eddington, and low black hole masses (e.g., 
Boroson & Green 1992, Pounds et al. 1995, Wang et al. 1996, Boiler et al. 1996, Laor et al. 
1997, Czerny et al. 1997, Marziani et al. 2001, Boroson 2002, Simkin & Roychowdhury 2003, 
Xu et al. 2003, Kawaguchi 2003, Wang & Netzer 2003, Grupe 2004, Grupe & Mathur 2004, 
Botte et al. 2004, Collin & Kawaguchi 2004). Apart from that, several other parameters 
were suggested to influence the main properties of NLSl galaxies, including orientation (e.g., 
Osterbrock & Pogge 1985, Bian & Zhao 2004), winds and density effects (e.g., Lawrence et 
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al. 1997, Wills et al. 2000, Bachev et al. 2004, Xu et al. 2003, 2006), metallicity (e.g. 

Mathur 2000, Komossa & Mathur 2001, Shemmer & Netzer 2002, Nagao et al. 2002, Warner 
et al. 2004, Fields ct al. 2005), and absorption (e.g., Komossa & Mccrschweinchen 2000, 
Gierlinski & Done 2004). There is evidence that NLSls follow a different black hole mass 
- velocity dispersion (Mbh — cr) relation than broad line Seyfert 1 galaxies (e.g., Mathur 
2001, Grupe & Mathur 2004, Mathur & Grupe 2005) with important implications for galaxy 
evolution and black hole growth. 

While the X-ray and optical properties of NLSls were explored intensively in the last 
decade, relatively little is known about their radio properties. Ulvestad et al. (1995) observed 
7 NLSl galaxies with the VLA and concluded that they are only of modest radio power and 
that radio emission is more compact than a few hundred parsec. Moran (2000) studied 24 
NLSls with the VLA and found that most of the sources are unresolved and show relatively 
steep spectra. Stepanian et al. (2003) found 9 radio(FIRST)-detected NLSls among 26 
NLSls of the Second Byurakan Survey, all of them radio-quiet. Greene et al. (2006) per- 
formed radio observations of 19 galaxies with low BH masses, several of them with optical 
spectra similar to NLSls. They only detected one source. 

So far, only a few radio- loud NLSls have been identified^. PKS0558-504 (Remillard et 

al. 1986, Siebert et al. 1999), RXJ0134-4258 (Grupe et al. 2000), SDSSJ094857.3+002225 
(Zhou et al. 2003) and SDSSJ172206.03+565451.6 (Komossa et al. 2006) all exceed radio 
indices R=10 and show optical NLSl spectra. RGB J0044+193 (Siebert et al. 1999) is 
radio-quiet most of the time (Maccarone et al. 2005), but may have been radio-loud at the 
epoch of the 87GB survey (Siebert et al. 1999, Maccarone et al. 2005) due to variability in 
the radio band. PKS 2004-447 (Oshlack et al. 2001) is a (non-typical) NLSl or possibly a 
narrow-line radio galaxy (Zhou et al. 2003) and is very radio-loud. 

It is still a puzzle why radio-loud NLSl galaxies are scarce. In particular, it is still 
unclear whether they are truly preferentially radio quiet, or rather the scarcity of radio- 
louds is due to some selection effect. This question can be addressed by a systematic search 
for, and study of, radio-loud NLSls. 

The study of radio properties of NLSls, in particular their radio loudness and radio 
variabihty, also allows us to re-address the question whether NLSl galaxies are preferentially 
viewed face on. Furthermore, radio observations of NLSl galaxies provide constraints on the 
couphng between jets and accretion disks (e.g., Zdziarski et al. 2003). There are indications 



^ there are preliminary reports that optical identifications of FIRST radio sources produce a higher rate 
of radio-loud NLSls than optical selection (Whalen et al. 2001; see our Section 5.8 for a comparison with 
Whalen et al. 2006) 
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that objects with accretion rates close to the Eddington rate tend to be radio-weaker; a 
behavior observed in Galactic X-ray binaries in soft/high-state and AGN (Maccarone et al. 
2003, Greene et al. 2006, and references therein). 

Here, we report results from a search for radio-loud and 'almost radio-loud' NLSl galax- 
ies, based on cross-correlating a catalogue of optically identified NLSls with radio catalogues. 
This is the first systematic study of (non-radio-selected) radio-loud NLSl galaxies, and of 
their radio, optical and X-ray properties. We use the term "Narrow-line Seyfert 1 (NLSl) 
galaxy" collectively for high-luminosity and low-luminosity objects, i.e. for narrow-line type 
1 quasars and for narrow- line Seyfert type 1 galaxies. A cosmology with Hq=70 km/s/Mpc, 
f2M=0.3 and Q\=0.7 is adopted throughout this paper. 

This paper is organized as follows: In Section 2 we describe our methods for cross- 
matching a large number of catalogues and present candidate radio-loud NLSl galaxies. In 
Sections 3 and 4 we provide results from our analysis of the optical and X-ray spectra of 
these NLSls. Consequences of our results for NLSl models, and mechanisms to explain the 
lower frequency of radio-louds among NLSls than quasars are discussed in Section 5. 

2. GAVO seeirch for radio-loud NLSl galeixies 

2.1. Search methods 

In order to search for radio-loud NLSl galaxies and study their properties, we extracted 
all NLSl galaxies included in the 11th edition of the "Catalogue of Quasars and Active 

NticIcI" compiled by Veron-Cetty & Veron (2003; henceforth referred to as VQC), separately 
for the catalogue of "faint" objects (i.e., Seyferts fainter than absolute magnitude M — 23) 
and the "bright" ones (i.e, quasars brighter than M = 23). We present results for the narrow- 
line type 1 quasars here, while results for the Seyfert galaxies will be reported elsewhere. The 
defining criterion for inclusion as NLSl galaxy in the VQC is a width of the H/5 emission 
fine FWHMh^ < 2000 km/s. The list of the 128 selected narrow-line quasars was then 
cross-correlated with the FIRST ( "Faint Images of the Radio Sky at Twenty centimeters" , 
at 1.4 GHz; White et al. 1997, Becker et al. 1995), NVSS ("NRAO VLA Sky Survey", at 
1.4 GHz; Condon et al. 1998), SUMSS ("Sydney University Molonglo Sky Survey", at 843 
MHz; Mauch et al. 2003), WENSS ("Westerbork Northern Sky Survey", at 0.33 GHz; de 
Bruyn et al. 1998), PMN ( "Parkes-MIT-NRAO" radio survey including all, the Southern 
Survey, Zenith Survey, Tropical Survey, and Equatorial Survey, at 4.85 GHz; Griffith et al. 
1994), 87GB ("Green Bank Radio Survey", at 4.85 GHz; Gregory & Condon 1991) and the 
PKS ("Parkes Radio Survey", at 2.7 GHz; Wright & Otrupcek 1990) radio catalogues, in 
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order to select all radio-detected NLSls, and among these the radio-loud ones. The input 
list was simultaneously cross-correlated with the USNO-Bl (Monet et al. 2003), USN0-A2 
and GSC2.2 optical photometry catalogues in order to obtain blue magnitudes. 

For cross- matching of radio catalogues with optical catalogues we used the cross- matcher 
application developed within the German Astrophysical Virtual Observatory (GAVO)^ project 
The GAVO cross-matcher is a multi-archive, multi-server, multi-catalogue, statistical spatial 
matcher operating solely on astrometric coordinates and their uncertainties (Adorf et al. 
2003, 2005). 

In our case, it takes as input an extract of the 'Veron quasar catalogue' which includes 
all NLSl galaxies in that catalogue and their coordinates. A query module then accesses 
the VizieR archive^ and extracts coordinates and flux and magnitude data from selected 
catalogues in the archive. 

While the VizieR match-list services are capable of deterministically matching the input 
list against a single catalogue, we then proceed in matching the different output lists against 
each other in order to find the most likely single multiwavelength "counterpart". In a 
first step, match-candidates are generated in the following way: combinations of potential 
counterparts are formed, which are viewed as hypothetical "clusters" on the sky. For each 
such cluster its central position is computed, followed by a computation of the projected 
distances on the sky between each counterpart candidate and the central cluster position. 

Next, based on the astrometric positions and associated uncertainties, for each coun- 
terpart in a given match candidate cluster its isotropic Mahalanobis distance (Mahalanobis 
1936, Devroye et al. 1996) to the central cluster position is computed. It measures the 
statistical distance to the cluster center. For each cluster a goodness-of-fit criterion is com- 
puted, namely the standard x^-metric, which is simply the sum of the squared Mahalanobis 
distances. It can be viewed as a measure of compactness of the cluster, where instead of the 
physical distances, the statistical distances are being used. 

Since the number of counterparts contributing to a match candidate generally differs 
from candidate to candidate, we cannot use the metric directly for comparing the goodness 
of the match. Instead the reduced metric is computed. For the required degrees of 
freedom A^f we use A^f = 2Ai — 2, where Ai denotes the number of counterparts contributing 
to a match candidate. We subtract 2 in order to allow for the estimation ("fitting") of the 
two sky-coordinates, i.e. right ascension and declination, when computing the cluster center. 



^Scc http://www.g-vo.org/ 

■^See http://vizier.u-strasbg.fr/viz-bin/VizieR 
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A reduced-x^ threshold parameter is used for discriminating against unreasonable match 
combinations. Finally, for each remaining match candidate, the cross-matcher extracts the 
astrometry and optical and radio photometry data. 

Position uncertainties, inherent to the different (radio) surveys, are either taken into 
account on an object-by-object basis if individual position errors are available (this is gen- 
erally the case for the FIRST and 87GB catalogues in VizieR; typically 1" (FIRST) and < 
30" (87GB)) or else typical average position imcertainties d for the different catalogues are 
used (explicitly, d=l" (USNO-A2,USNO-Bl,GSC2.2), 7" (NVSS), 30" (PKS and PMN), 10" 
(SUMSS and WENSS)). The catalogue-specific search radii were typically larger by a factor 
of 3. 



2.2. Safety checks and data screening 

While this automated way of simultaneously cross-correlating various catalogues and of 
selecting counterparts is generally very successful, there are still a number of "peculiarities" 
inherent to the way the catalogues were produced in general, or to individual galaxies in 
particular, which affect the cross-matching output and which make it necessary to check a 
number of results "by hand" . While the long term goal is to include most of these effects 
into the cross-correlation software, in this pilot project we paid special attention to actually 
identifying potential problems and effects that need closer attention in the future when 
applying similar cross-correlation methods to much larger samples. 

Nearby bright galaxies are missed in blue magnitude catalogues because they are too 
extended and too bright for standard automated magnitude measurements to still work. 
This is not a problem for the present sample since the quasars are sufficiently distant but 
will be more severe for samples of nearby Seyfert galaxies. 

A few catalogued NLSls still have optical coordinate off-sets which are larger than 
expected in the optical band. One reason is that their IRAS coordinates rather than optical 
coordinates were still listed, which come with a larger uncertainty. Therefore, counterparts 
in blue magnitudes and FIRST radio sources might be missed. For the present sample, 
IRAS11598-0112 and RXJ01354-0426 were not assigned optical counterparts by the original 
cross-correlation procedure. IRASl 1598-01 12 had its IRAS coordinates catalogued, while 
the true optical coordinate then turns out to be consistent with the FIRST radio position. 
The optical counterpart of RXJ01354-0426 is off-set by ~3-4" from the coordinate in the 
VQG. However, since there is no other good multi-wavelength counterpart in the field, we 
assume that optical, radio and X-ray source are one and the same. 
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Bright radio sources are often detected in several surveys with widely different positional 
uncertainties. In that case it is generally a good idea to then use the smallest positional error 
among the radio catalogues for further cross-correlation with non-radio catalogues to reduce 
the number of potential mult i- wavelength match candidates. In the present study, in two 
cases, a bright source was detected in several radio surveys, suggesting wc deal with one and 
the same source, but the source coordinates in those radio surveys with small coordinate 
uncertainties were no longer consistent with the optical position. In the specific case of 
RXJ07101-I-5002, only a 87GB counterpart was found by cross-correlation which is 45" off 
the optical position. The 87GB source has counterparts in several other radio catalogues 
including NVSS and FIRST, which show the same off-set from the optical coordinates and 
similar radio fluxes. This means the radio source is one and the same in all radio catalogues. 
However, given the small coordinate uncertainty in the FIRST catalogue, it can no longer be 
related to the optical NLSl galaxy. The FIRST image reveals a complex radio structure with 
several bright knots and lobes. We have checked that none of these sub-structures coincides 
with the optical NLSl. Rather, the core of the radio galaxy does have a 'normal' galaxy 
as optical counterpart, reported by Bauer et al. (2000). The question remains whether the 
X-ray source RXJ07101-I-5002 then actually is the counterpart to the radio source (as often 
assumed; e.g., Bauer et al. 2000), or whether it is the counterpart to the optical NLSl galaxy, 
as concluded by Xu et al. (2003) who provided the optical spectroscopy of that source. The 
radio galaxy is radio loud with i?i.4=64. We exclude it from the sample, though, since the 
NLSl galaxy is not its counterpart. In the case of SBS1330+519, only a WENSS counterpart 
was found, off-set by 20". The WENSS source has a FIRST and NVSS counterpart which, 
however, are too far off from the optical position of SBS1330-I-519 to be still considered as 
counterparts. SBS1330-I-519 was thus excluded from our sample. We note in passing that 
the relatively strong radio source itself does have a very faint optical counterpart in the 
Sloan Digital Sky Survey (SDSS; York et al. 2000) photometric catalogue (Abazajian et al. 
2005), SDSSJ133219.38+514439.9. With SDSS r and g magnitudes of 21.5 and 22.3 mag, 
respectively, this background object is very radio loud (i?i.4=8940). No optical spectrum is 
available. 

Radio images of all radio-detected NLSls were inspected by eye to convince ourselves 
of the reliability of counterpart identification. 

2.3. Assessment of radio loudness 

For all confirmed radio-detected objects their radio loudness was then estimated, fol- 
lowing Kellermann et al. (1989), who define the radio index R as ratio of 6cm radio flux to 
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optical flux at 4400A. i? = 10 is commonly used to mark the "border" between radio-quiet 
and radio-loud objects, originally calculated based upon the assumption of similar spectral 
shapes in the optical and radio band with a = —0.5 (Kellermann ct al. 1989). Radio- 
surveys employed in the present study generally measure at 1.4 GHz. Still keeping the same 
assumptions, in particular, a = —0.5, we compute the radio index at 1.4 GHz as -Ri.4=1.9i?, 
and we use i?i.4=19 as a reasonable distinction between radio-loud and radio-quiet objects. 
Optical blue magnitudes mB2 were extracted from the USNO-Bl catalogue and used to cal- 
culate RiA- For comparison, blue magnitudes from the first epoch of the USNO-Bl (mei), 
the GSC2.2 (mej) and the USN0-A2 (me) catalogues were also collected to get an impres- 
sion on the uncertainty in blue magnitude, which would be either due to source-intrinsic 
variability and/or measurement uncertainties. 

As a last step, only a careful object-by-object check then allows to scrutinize which 
objects are actually radio-loud NLSl galaxies, or good candidates. In particular, peculiar 
intrinsic spectral shapes, high redshift, or heavy absorption, which will affect radio and 
optical band in a different way, would invalidate the simple way of calculating the radio 
index R as ratio of observed radio to optical flux density. Also, merging galaxies with 
close double nuclei make counterpart identiflcation ambiguous. Finally, optical spectroscopic 
classiflcation of each radio-loud candidate needs to be checked carefully in order to conflrm 
its NLSl nature. An uncertainty which generally cannot be overcome with existing data 
given their non-simultaneity, is source variability either in the optical and/or radio band. 



2.4. Sesirch results 

2.4-1- General results 

Among the 128 NLSls in the VQC we flnd that 35 have radio counterparts in one or 
more catalogues. In Fig. 1 we present the distribution of their radio indices Ria- Formally, 
among this sample, ~11 objects exceed or are close to the "critical" value of i?i.4=19 com- 
monly used to distinguish between radio-quiet and radio-loud objects. The radio-loud NLSls 
are TEXllllH-329, SDSSJ094857.3+002225, SDSSJ172206.03+565451.6, 1RAS20181-2244, 
RXJ0134-4258, PKS0558-504, 1RAS09426+1929, SBS1517+520, RXJ16290+4007, 2E1346+2637, 
and RXJ23149-I-2243 (Tab. 1; sources are listed by their name as given in the VQC). Among 
these, we recover the previously known radio-loud NLSl galaxies PKS0558-504, RXJ0134- 
4258, SDSSJ172206.03+565451.6 and SDSS J094857.3+002225.^ PK2004-447, not included 



^Wang et al. (2004) claimed the detection of another radio-loud NLSl, SDSSJ022119.8-F005628.4. How- 
ever, the position off-set between SDSS galaxy and the faint NVSS radio source is 24", much larger than the 
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in our sample since it is a Seyfert rather than a quasar, is unusual in almost lacking Fell 
emission in its optical spectrum (Oshlack et al. 2001), which led to the suggestion it could be 
a narrow-line radio galaxy or a type-2 AGN. Given the large scatter in Fell of the radio-loud 
objects of our sample, it is well possible that PK2004-447 is similar, just at the low Fell end 
and we include it in Tab. 1. 

Based on average blue magnitudes of the non-detected sources and given the NVSS flux 
limit (~2.5mJy), no more very radio-loud sources are expected among the non-dctcctions. 
Several more borderline objects close to -R1.4 ~ 19 could exist, but only if their radio-flux 
was just exactly below the NVSS flux limit. 

Radio source positions agree well with the optical coordinates of the sources. FIRST 
counterpart position off-sets arc typically <1", while NVSS off-sets arc typically <5". The 
radio emission of the radio-loud NLSl galaxies is generally compact. Their FIRST extent is 
always <2" which is consistent with point sources (White et al. 1997). 

The 35 radio-detected sources cover about four orders of magnitude in radio index R. 
Below log-R ~ 1 (Fig. 1) the data hint at a rather smooth distribution of radio indices, 
while only few objects are well above log/? ~ 1. A number of our objects (Fig. 2) fall into 
that part of the radio-loud regime sometimes referred to as 'radio-intermediate' (logi? ~ 
1-2, as opposed to 'truly radio- loud' above logi? ~ 2 — 3; e.g., Falcke et al. 1996), a regime 
where objects are relatively rare, and which classically defined the radio-loud radio-quiet 
bimodality of AGN. 

2.4-2. Variability 

Several sources were observed during both the NVSS and FIRST survey, and these 
data can be used to assess the radio variabihty of the galaxies. Generally, we find that 
NVSS and FIRST radio fluxes agree well with each other. For three sources, fluxes differ 
by approximately a factor of 1.6-1.8. Among these, variability of SDSSJ094857.3+002225 
(integrated flux density -S'first = 111.5 mJy and S'nvss = 69.5 mJy) is likely real (Zhou et 
al. 2003). The other two sources are Mrk 478 (5'first = 3.3 mJy and S'nvss = 5.2 mJy) 
and PHL1811 (^first = 1-2 mJy and ^nvss = 2.1 mJy). Higher NVSS than FIRST fluxes 
are actually expected in case of very extended sources because part of the extended emission 



expected off-set of ~7" for faint sources. Furthermore, wc do not find any radio counterpart in tfie FIRST 
survey, indicating that the very faint NVSS counterpart, if real, is either very extended or variable. Deeper 
radio observations are needed to confirm the existence of this radio source, and to see whether it indeed is 
the counterpart to the NLSl galaxy. 
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will be missed by FIRST (e.g., Fig. 7 of White et al. 1997). Indeed, we find indications 
that the radio emission of PHL1811 is extended (major axis of 2.9" reported in the FIRST 
catalogue, and a value of < 74" in the NVSS catalogue). Inspecting the NVSS and FIRST 
images wc find that Mrk 478 has a nearby bright companion source which might have slightly 
contributed to the NVSS radio flux estimate of Mrk 478. Barvainis et al. (1996) report a 
dedicated radio observation of Mrk 478 at 1.49 GHz with S = 4.46 ± 0.47 mJy, within the 
errors almost consistent with the NVSS and FIRST measurements. 

Among those sources detected in the radio band there are 10 in the FIRST data base 
which do not have a NVSS counterpart, while 7 in the NVSS data base do not have a 
FIRST counterpart. The latter is just due to the non- availability of FIRST observations at 
the respective coordinates. The radio images of sources with FIRST-only detections were 
carefully inspected by eye and it turns out that all FIRST detections are consistent within 
better than a factor of two with the NVSS upper hmits. 

2.4-3. Radio spectral indices 

Several sources were observed in the radio regime at more than one frequency. These 
data can be used to estimate radio spectral shapes. Radio spectral indices are reported in 
Tab. 1. shows a large scatter. Steep and flat spectrum sources are observed. While the 
two galaxies SDSSJ094857.3+002225 and RXJ16290+4007 show very flat spectra (0:^=0.6 
and 0.4, respectively) the majority of NLSl galaxies are of steep spectrum type (ctr < —0.5). 

2.4.4- Notes on individual objects 

We comment here on the 11 apparently radio-loudest objects of this study, i.e. those 
faUing above or close to the critical radio index, Ria — 19 (Fig. 1). Emphasis is put on the 
question how safely these galaxies can be called "radio-loud" (Sect. 2.3) and how rehable 
their optical spectral classiflcation is. 

TEXl 1111+329 is an ultraluminous infrared galaxy (ULIRG). Optical spectroscopy 
was provided by Zheng et al. (2002) who classify TEXlllll+329 as NLSl galaxy (FWHMh/? ~ 
1980 km/s, [OIII]/H/3=0.75, FeII/H/3=l.l). Lipari et al. (2003) report the detection of a 
high- velocity outflow component in [OIII]5007 with v ~ 1300 km/s. The galaxy is heavily 
absorbed, with at least ii^B-v=l-08, estimated from the observed Balmer-line ratio (Zheng 
et al. 2002). Assuming that this obscuration is not intrinsic to the line-emitting clouds 
(the BLRs in NLSls might be sufficiently distant from the nucleus that dust can survive) 
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and correcting the observed optical continuum spectrum by Ab=4.5 leads to a much lower 
value of Ri 4 ~ 20. On the other hand, TEXlllll+329 would also be classified as radio 
loud just based on its radio power, which is Pi 4 ~ 10^^ W/Hz. The radio spectrum declines 
toward higher frequencies, and TEXlllll+329 is not detected in an imaging survey at 15 
GHz (Nagar et al. 2003). We keep the source as candidate radio-loud NLSl. 

SDSSJ094857.3+002225 is a very radio-loud NLSl galaxy (Zhou et al. 2003), vari- 
able in the radio and likely the optical band. Inspecting the NVSS and FIRST radio image, 
we noticed a second radio source very close to SDSSJ094857.3+002225, F1RSTJ094901. 5+002258, 
at a separation of 71" (corresponding to a projected distance of 395 kpc) and extended 
by ~5". Suspecting it to be either a knot in a radio jet of SDSSJ094857. 3+002225, or a 
background source, we had a closer look at optical images and cannot identify any optical 
counterpart. Closest is a faint SDSS galaxy, SDSSJ094901.9+002306.4, at an off-set of 10", 
an unhkely counterpart of the FIRST radio source. It is possible that the radio source has 
a faint, so far undetected, optical counterpart. 

SDSSJ172206.03-I-565451.6 is a very radio-loud NLSl galaxy. Depending on choice 
of blue magnitude, it may be comparably radio-loud as SDSS J094857. 3+002225. Its radio 
loudness and optical spectral properties are discussed in detail by Komossa et al. (2006). In 

particular, its black hole mass in unusually small given its radio loudness (Komossa et al. 
2006). SDSSJ172206. 03+565451. 6 is the optically most variable among our sample. 

RXJ0134-4258 is a radio-loud NLSl galaxy (Grupe et al. 2000) with pecuhar X-ray 
spectral variabihty (Komossa &: Meerschweinchen 2000, Grupe et al. 2000). It was detected 
in the PMN survey and its radio-loudness was confirmed by a dedicated follow-up observation 
by Grupe et al. 

IRAS 20181-2244 was identified as a quasar with narrow emission fines by Elizalde 

&; Steiner (1994). Halpern & Moran (1998) rejected a "type-2 quasar" classification of 
IRAS20181-2244 and called it a NLSl-type object, with FWHMh/3=460 km/s and the pres- 
ence of broad wings in the Balmer lines. Kay et al. (1999) found low-level continuum 
polarization (2%) but no broad Balmer lines. The X-ray emission of IRAS20181-2244 (e.g., 
Vaughan et al. 1999) is rapidly variable (Halpern & Moran 1998, Leighly 1999). The Balmer 
decrement of the narrow lines implies significant reddening (Tab. 1 of Halpern & Moran 1998) 
but the rapid X-ray variability argues against a full obscuration of the continuum source it- 
self. The ASCA X-ray spectral fits indicate some excess absorption along the line of sight 
and Halpern Sz Moran give £'b-v ~ 0.3 mag as a reasonable estimate of extinction intrinsic 
to IRAS 20181-2244. In the radio band, IRAS 20181-2244 is also detected during the WISH 
(Westerbork in the Southern Hemisphere) survey at 352 MHz (de Breuck et al. 2002). 
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RXJ16290+4007 is a bright source in the X-ray (e.g., Giommi et al. 1991, Elvis et 
al. 1992, Brinkmann & Siebert 1994) and radio band (e.g., Giommi et al. 1991, Laurent- 
Muchlciscn ct al. 1997, White et al. 2000) likely variable by ~0.5 mag optically (Helfand 
et al. 1991). RXJ16290+4007 is known as fiat spectrum radio quasar (FSRQ). Its FSRQ 
nature, with respect to its X-ray and radio properties, is discussed in detail by Padovani et 
al. (2002). Steady TeV gamma ray emission from RXJ 16290+4007 was searched for, but 
not detected, with the Whipple telescope (Falcone et al. 2004). Optical spectra provided by 
Bade et al. (1995) hinted at the NLSl nature of RXJ16290+4007, which was confirmed by 
Schwope et al. (2000) and Grupe et al. (2004). 

While Padovani ct al. (2002) concentrated on the radio and X-ray properties of RXJ 16290-1-4007, 
interpreted its relatively steep X-ray spectrum as due to synchrotron radiation, and sug- 
gested RXJ16290-I-4007 is the first member of a newly established class of HBL(high energy 
peaked BL Lac)-like FSRQs, we rather follow the alternative idea that the optical and X-ray 
properties of this source are driven by its "narrow- line Seyfert 1 character" . 

IRAS09426-I-1929 is ultraluminous in the infrared with a single nucleus. Zheng et 
al. (1999) identify as counterpart of the IR source a NLSl galaxy at 2;=0.284, based on 
an optical spectrum which shows weak [OIII] , strong Fell and narrow Balmer lines (no hue 
parameters were measured). On the one hand, the host galaxy is detected as well and 
may contribute to the blue magnitude which would make the computed value of the radio 
index R a lower limit. On the other hand, extinction measurements arc not yet available 
for 1RAS09426+1929. Since the cores of ultraluminous infrared galaxies arc often heavily 
obscured, any calculation of R based on the observed optical spectrum would then make the 
value of R an upper limit. This second effect is hkely to exceed the first. 

PKS0558-504 is a known radio-loud source (Rcmillard et al. 1986, Siebert et al. 1999), 
with an optical NLSl spectrum (Rcmillard ct al. 1986). It is highly variable in the X-ray 
band (e.g., Rcmillard et al. 1991, Gliozzi ct al. 2001, Wang et al. 2001). The rapid X-ray 
variabihty and the lack of excess X-ray absorption (O'Brien et al. 2001, Brinkmann et al. 
2004) indicate an unobscured view on the central engine. 

RXJ23149-|-2243 was identified as quasar with strong Fell emission by Wei et al. 
(1998). Xu et al. (2003) report emission-hne measurements which indicate a relatively 
large width of the Balmer lines. Our newly taken spectrum (Sect. 3) puts the galaxy at the 
border between NLSl and Seyfert 1, depending on the way the width of H/3 is measured. The 
[OIII] 5007 hue of this galaxy shows an exceptionally strong blue wing which is blueshifted 
from the narrow core by 1260 km/s (Sect. 3). RXJ23149-F2243 is detected by IRAS and 
turns out to be a luminous infrared galaxy (LIRG). 



-14- 



2E1346+2637 is an ultrasoft AGN (photon index ~ -3.7) at high redshift, z=0.92 
(Puchnarewicz et al. 1994), optically classified as NLSl galaxy (FWHMh/3 — 1900 km/s, 
weak [OIII] and presence of Fell; Puchnarewicz et al. 1992). The X-ray spectrum of the 
galaxy displays one of the strongest soft excesses known. Puchnarewicz et al. (1994) report 
a radio upper limit of 6 mJy at 4.85 GHz and mention that this leaves a small margin that 
2E1346H-2637 is radio-loud. 2E1346+2637 is detected by FIRST but not NVSS. The NVSS 
upper limit is consistent with the FIRST measurement. 

SBS1517-|-520 is listed as NLSl in the VQC but little else is known about this galaxy 
(Bicay ct al. 2000 reported its optical position). Our analysis of the SDSS spectrum (Sect. 
3) does not confirm its classification as a NLSl galaxy in the sense that FWHMH/jg = 
3220 km/s (based on single-component Gauss fit (Tab. 2); only the 'direct' determination 
of the width of FWHMh/Jj = 2030 km/s, puts it close to the border). Because other 
emission-line parameters share similarity with those of NLSl galaxies, we continue to include 
SBS1517+520 in our tables. 



3. Optical spectroscopy of the radio-loud NLSls 

In order to measure the optical emission-line and continuum properties of the radio- 
loud NLSls in a homogeneous way, we (re-)analyzed their optical spectra when available. 
The spectrum of RXJ0134-4258 was taken by one of us (D.G.) at ESO's 1.5m telescope 
(Grupe 2004), as was the spectrum of RXJ23149+2243 (D.X.), at the Xinglong 2.16m tele- 
scope (Xu et al. 2003), while those of SDSSJ0948573-f002225, SDSSJ172206.03+565451.6, 
RXJ16290+4007, SBS1517+520 were extracted from the SDSS data basc^ In addition, a 
new spectrum of RXJ23149+2243 was secured on November 25, 2005 at the 2.16m Xinglong 
telescope, motivated by the remarkable structure of the [OIII]5007 profile of this galaxy. 
SDSS spectra of the other galaxies are not (yet) available. For them, emission-hne data were 
collected from the literature when pubhshed (Tab. 2). 

The SDSS DR4 products were used. These spectra arc flux- and wavelength-calibrated 
by the spectroscopic pipeline in the course of the DR4 processing. The spectrum newly 
taken with the 600gmm~^ grating at the 2.16m Xinglong telescope was reduced following 
standard procedures (see Xu et al. 2003). The IRAF package SPECFIT (Kriss 1994) was 
used for spectral analysis. 

We first corrected the spectra for Galactic extinction, using the Galactic Eb-v in the 



^See http://cas.sdss.org/dr4/en/tools/explore/obj.asp 
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direction of the individual galaxies (Schlegel, Finkbeiner & Davis 1998) and a R=3.1 ex- 
tinction law. We then fit a power law to the underlying continuum using the "continuum 
windows" at 3010-3040, 3240-3270, 3790-3810, 4200-4230, 5080-5100, 5600-5630, 5970-6000, 
and 6005-6035 A (Forster et al. 2001, Vanden Berk et al. 2001) known to be relatively 
free from strong emission lines. The continua are generally well fit by single powerlaws with 
indices between a = —1.3 and —2.3 (/a oc A"*""). 

Since the optical-UV spectra of all sources show emission from Fell complexes, an 
Fell spectrum was then fitted and subtracted (excluding the UV Fell multiplets around 
MgIIA2798 which are not included in the template) using the optical Fell template and the 
technique described by Boroson & Green (1992). 

The Fell- and continuum-subtracted spectra were then used to measure emission-line 
properties. The emission lines were modeled by single Gaussian profiles or combinations 
of Gaussians and Lorentzians. The approach was to fit single-component Gaussians to the 
weak forbidden lines, apply two-component Gaussians to [OIII] when necessary, and compare 
double Gauss and Gauss-plus-Lorentz profiles in their success in representing the Ha and H/3 
line profile. In addition, a direct measurement of the H/3 line width was performed (without 
any assumption on line profile shape). H/3 line results are marked with an index 'g' or 'd' if 
referring to a Gaussian fit or the 'direct' method, respectively, and with 'b', 'n' or 'totl' when 
referring to the broad, narrow or whole emission component, respectively, of a two-Gaussian 
decomposition of the line profile. 

Results for the diagnostically most important emission fines are fisted in Tab. 2. For 

the sake of simplicity and homogeneity, the listed emission-line ratios are based on single 
Gaussian fits to the lines, even though for some sources this is not a perfect match of the 
line profile. The instrumental resolution was corrected for in the FWHMs measured by 
us. We confirm the NLSl classification for almost all galaxies, with the possible exceptions 
of SBS1517+520 and RXJ23149+2243 (see below). Line widths of confirmed NLSls range 
between FWHMn/jg = 930-1580 km/s, while the fiux ratio [0III]/H/3g varies between 0.04 
-1.1 (results from Gauss fits are reported here). Fell is detected in all spectra and ranges 
between FeII4570/H/3g = 0.5-3.2 (FeII4570/H/3toti = 0.5-2.4). The NLSl galaxies of our 
sample cover almost the whole range of FeII4570/H/3 observed in the NLSl population (e.g.. 
Tab. 3 of Veron-Cetty et al. 2001), except that the lowest FeII/H/3 ratios are missing in our 
sample. 

SBS1517-F520 and RXJ23149-F2243 have rather broad Balmer fines, and it depends on 

the way the profile of H/3 is fit whether or not they formally fullfil the criterion of FWHMh/j < 
2000 km/s. Both of them have FWRMh/j^ > 2000 km/s while FWRMh/j^ < 2000 km/s. 
While the strict cut in FWHM is a useful working criterion for extraction of NLSl galaxy 
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candidates from larger samples, it has been repeatedly noted that the cut value is not well 
defined or even completely arbitrary. Veron-Cetty et al. (2001) suggested that a better 
classification criterion than FWHMh/3 may actually be the Fell ratio, FeII4570/H/9toti > 0.5. 
According to this criterion, both SBS1517+520 and RXJ23149+2243 still classify as NLSls. 
After these cautious comments, wc continue to keep the two galaxies in the sample, even 
though they fail a strict FWEMh/j^ < 2000 km/s criterion. 

RXJ23 149+2243 is exceptional in showing a very strong and broad blue wing in the 
[OII1]5007 and [0111]4959 lines. The presence of the blue wing was hinted in the earlier 
low-resolution spectrum of Xu ct al. (2003) and was confirmed in the new spectrum we 
took in November 2005 (full results on this galaxy will be presented elsewhere). [OIII]5007 
is well modeled by two Gaussian fines. The core, with FWHM[oiii]„=620 km/s, is at the 
same redshift as the peak of H/3 while the wing, with FWHM[oiii]|,=1560 km/s, is highly 
blueshifted by 1260 km/s. SBS1517+520 also shows a blue wing in [OIII]. The fine is weU 
fit by two Gaussians with FWHM[oiii]„=610 km/s and FWHM[oiii]b=2070km/s, blueshifted 
by v=90 km/s. In the SDSS spectra of the other objects of our sample the blue wing is 
present, but not exceptionally strong. For SDSSJ172206. 03+565451. 6 the line blueshift is 
v=UO km/s, for RXJ16290+4007 we measure ^=230 km/s. 



4. X-ray spectroscopy and veiriability 

While some of the sources of our sample have no previous X-ray studies, others (in 
particular, PKS0558-504) have been observed by most major X-ray satellites. In order to 
obtain representative X-ray fluxes of all our sources in a homogeneous way, including upper 
limits for non-detected sources, we make use of the data obtained with the PSPC (Positional 
Sensitive Proportional Counter) during the ROSAT all-sky survey (RASS; Voges et al. 1999) 
and the later phase of pointed observations. 

9 out of the 11 soures have X-ray detections, 8 of these in the RASS and 3 in pointed 
observations. Countratcs and hardness ratios were obtained for all sources. Fixing the ab- 
sorption at the Galactic values in the directions of the individual sources, hardness ratios 
were then converted to powerlaw indices, and the absorption-corrected fiux in the (0.1-2.4) 
keV band was determined. Powerlaw indices are in the range Fx = —1.9.. —4.3 (Tab. 2) with 
the exception of the faint source SBS1517+520 (Fx = —1-1, fikely caused by excess absorp- 
tion). In case pointed observations existed (RXJ16290+4007, 2E1346+2637, RXJ0134-4258) 
these were given preference over the RASS data. 

For the two non-detected sources, TEXlllll+329 and IRAS09426+1929, upper limits 
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were obtained. However, we caution while soft X-ray non-detections convert into safe upper 
limits for non-obscured sources, the same is not true for obscured sources. This is likely the 
case for both TEXll 11 1+329 and IRAS09426+1929. For completeness, we report the 2a 
upper limits on the soft (0.1-2.4 keV) X-ray count rates in Tab. 2, but do not convert them 
into luminosities, and do not use them for accretion rate estimates. 

We also inspected the RASS X-ray lightcurves of the X-ray brighter sources and find 
that most are consistent with constant source flux except PKS0558-504 which shows repeated 
flaring activity. Variability of the X-ray fainter sources cannot be judged due to the poor 
statistics. 



5. Discussion 

Judged on radio index, 11 objects classify as radio-loud NLSl galaxies or candidates. 
Our study almost triples the number of recognized radio-loud NLSls if the candidates close 
to the deviding line are confirmed. Implications of the new results are discussed below. 

5.1. Compcirison of the properties of the known radio-loud NLSl galeixies 

Even though the number of known radio-loud NLSls and the candidates is still small, 
we have enough objects now to compare their X-ray, optical and radio-properties, to have a 
first assessment of their similarities or differences. In particular, we checked whether the ~11 
radio- loud NLSls are extreme in one or more of their multi- wavelength spectral properties, 
i.e., whether they differ significantly from the radio-quiets in their properties, and whether 
there is a trend within the sample such that one or more properties correlate with radio 
loudness. 

Among optical emission lines, the "first Eigenvector" (Boroson & Green 1992) is evident 
among the objects of our sample in that weak [OIII]/H/? emission tends to go hand in hand 
with strong Fell emission and vice versa. Both, [OIII]/H/? and Fell//? vary strongly across 
the sample, almost spanning the whole range observed in NLSl galaxies. In Fig. 3 we plot 
the location of the radio-loud NLSls in the FeII-[OIII] diagram, in comparison with one of 
the largest homogeneously analyzed optical NLSl samples available to date (Xu et al. 2006). 
We do not find a strong correlation of any of the parameters with radio loudness, but larger 
samples are needed to confirm this. 

Even though radio loudness did not occur in the sample of Marziani et al. (2001) for fine 
widths FWHMh/3i, < 4000 km/s, the radio- louds of our sample do populate similar regions 
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as other AGN in the FWHMH^j^-Fell diagram (Fig. 4 of Marziani et al. 2001), and extend 
the range where radio loudness is observed to significantly smaller FWHMs than previously 
known (see also Sect. 5.5). 

Given suggestions that the blue wings in the [OIII]5007 hne of quasars are sometimes 
hnked to interactions of the radio plasma with the NLR clouds (Leipski & Bennert 2006, 
and references therein) we paid special attention to the existence and strengths of blue wings 
in dependence of radio loudness but do not find any correlation, even though 2-3 out of the 
11 sources do show very intense blue wings, two with strong blueshifts, among the highest 
values ever observed. 

Most NLSls of our sample have steep radio spectra (Tab. 1) in the frequency band 0.33- 
1.4-4.85 GHz^ and their radio emission is compact. As such, they share some similarity with 
'compact steep spectrum sources' (CSS, Peacock & Wall 1982) which are believed to be 
young radio sources which did not yet evolve strongly (Marecki et al. 2005, and references 
therein). SDSSJ094857.3+002225 and RXJ16290+4007 have inverted radio spectra. Only 
one source, SDSSJ094857.3+002225 (Zhou et al. 2003), shows significant variability in the 
radio-band by a factor ~1.6. 

Soft X-ray spectra are generally steep, but come with large errors. The mean ROSAT 
powerlaw index of the radio-loud NLSls of our sample'' (excluding SBS1517-I-520 which has 
unusually flat index). Fx = —2.7, is similar to that of soft X-ray selected NLSl galaxies 
(Fx = -2.96 ± 0.41, Grupe et al. 2004). PKS0558-504 (e.g., Remillard et al. 1991), 
IRAS20181-2244 (Halpern & Moran 1998) and RXJ16290+4007 (Padovani et al. 2002) are 
highly variable in the X-ray band while RXJ0134-4258 strongly changed its spectral shape 
(Grupe et al. 2000, Komossa & Meerschweinchen 2000). 

5.2. Starburst contribution ? 

It was suggested that NLSls are AGN in an early phase of evolution. This includes 
the possibility of enhanced starburst activity in these objects (e.g., Mathur 2000, Shemmer 
et al. 2004). Could starbursts contribute to or dominate the radio power of the radio- loud 
NLSls of our sample? Generally, normal and starburst galaxies are much less luminous than 



®Thc non-simultaneity of the data, taken with different instrumental configurations, has to be kept in 
mind, but unlikely all scatter in among the sample can be traced back to variability, else variability in 
individual cases would be huge. 

''note that a number of them were first identified by X-ray observations 



-19- 



quasars in the radio regime. In particular, all of our NLSls also exceed the radio powers of 
a sample of the most radio-luminous starbursts, all of them in luminous and ultraluminous 
infrared galaxies (Smith et al. 1998), logP4 85 SB — 22.3 — 23.4 W/Hz. Radio powers place 
all of our NLSls with i?i 4 > 19 in the radio-loud regime (logP > 24, Joly 1991), some of 
them close to the lower border, though. 

At first glance, a significant starburst contribution might be suspected for the two ultra- 
luminous infrared galaxies among our sample, TEXlllll+329 and 1RAS09426+1929. We 
used the tight correlation between IR luminosity and radio luminosity of starburst galaxies 
(Yun et al. 2001, and references therein) to predict the expected starburst radio power, 
given the 60/x IR luminosities of these two galaxies measured by IRAS. According to equ. 
(4) of Yun et al. (2001), with logPi,4,Nvss - 25 W/Hz TEXlllll+329 exceeds by a factor 
of 10 the contribution expected from a starburst, while for IRAS09426-I-1929 the measured 
radio power, logPi.4,Nvss — 24.7 W/Hz, is a factor 5 above the predicted starburst contri- 
bution. TEXll 11 1-1-329 has radio observations at several frequencies. The spectral index is 
relatively steep, but comparable to other sources of our sample. The spectral index between 
0.33-1.4 GHz is similar to SDSSJ172206.03+565451.6 and PKS0558-504 while towards higher 
frequencies it is comparable to RXJ0134-4258. 

The same comparison between observed radio power and expected starburst contribution 
was done for the three other galaxies with IRAS 60fi flux measurements. We find that the 
observed radio powers of PKS0558-504, IRAS20181-2244, and RXJ23149+2243 exceed by 
factors of 115, 6, and 6, respectively, the predicted starburst contribution. 

We conclude that radio-luminous starbursts are unlikely to dominate the radio emission 
of most of our sources. 



5.3. Black hole masses 

Reverberation mapping of Seyfert galaxies established a relation between the radius of 
the broad line region and the optical luminosity (Kaspi et al. 2005, Peterson et al. 2004, and 
references therein). Assuming that the broad line region (BLR) clouds are virialized (e.g., 
Wandel et al. 1999), the black hole mass can then be estimated as Mbh = G^^ Rblrv"^- 
The BLR radius, Rblr, in dependence of the optical luminosity at 5100A is given by 

Pblr = 22.3[^^gg^]0-69 Id (equ. (2) of Kaspi et al. 2005). The velocity v of the BLR 

clouds is usually estimated from the FWHM as v = f FWHM where / = ^ for an isotropic 
cloud distribution. While the relation is now well established for nearby Seyfert galaxies, few 
NLSls have been reverberation mapped so far (Peterson et al. 2000). Following standard 
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practice, we assume equ. (2) of Kaspi et al. (2005) is indeed applicable to NLSl galaxies as 
well and proceed in estimating the BH masses of our sample. 

La(5100A) was derived from the SDSS spectra (corrected for Galactic extinction) when 
available, or else estimated from blue magnitudes mB2 of the USNO-Bl catalogue (note that 
some of the galaxies of our sample vary in the optical band which introduces an uncertainty 
in the black hole mass estimates. Most objects, however, vary by less than a factor 2). The 
FWHM of the broad component from the two-component fit of H/9 was used for the black 
hole mass estimates reported in Tab. 3 when optical spectra were available. Else, results on 
FWHM from single-component line fits collected from the literature (Tab. 2) were used for 
an ordcr-of-magnitudc estimate. 

We obtain black hole masses in the range Mbh — (2 — 10) 10^ Mq (objects with two- 
component decomposition of H/5 available) and Mbh — (0.2 — 9) 10^ Mq (objects with single 
component fits to H/3; this underestimates FWHMH/3b and therefore underestimates the black 
hole mass) for our sample. These masses are surprisingly small, given the radio loudness 
of the galaxies (Fig. 4). In fact, all galaxies fall right into an unpopulated regime of the 
"Laor diagram" (Fig. 2 of Laor 2000) which plots radio loudness in dependence of black 
hole mass. While larger samples filled up some originally empty areas in the Laor diagram, 
and while some discussion has emerged whether there is any clear dependence of black hole 
mass on radio loudness at all (e.g., Lacy et al. 2001, Oshlack et al. 2002, Woo and Urry 
2002, McLurc & Jarvis 2004, Best et al. 2005), our radio-loud NLSl galaxies arc still located 
in an area very sparsely populated by galaxies. Firstly, because the black hole masses are 
unusually small given the radio loudness of the galaxies of our sample (e.g.. Fig. 2 of Laor 
et al. 2000, Fig. 3 of Lacy et al. 2001, Fig. 2 of McLure & Jarvis 2004, Fig. 10 of Metcalf 
& Magliocchetti 2006). Secondly, because most radio indices R of our objects are in a range 
which is classically known to be less populated by quasars and which defined the radio-loud 
radio-quiet bimodality of AGN (with radio-quiet quasars below it! ~ 1 — 10 and many radio- 
loud quasar above R ^ 100 — 1000 (e.g. Fig. 10 of Woo & Urry 2002)). An even more 
extreme outlier is PKS2004-447 which shows lower BH mass and higher radio index R than 
objects of our sample (Oshlack et al. 2001). 

While of unusually low mass when compared to radio-loud quasars, among the NLSl 
population itself out radio-loud NLSls are actually at the high-mass end (for comparison, in 
the sample of Grupe & Mathur (2004) NLSl masses range between 6 10^-4 10^ Mq with only 
one object around 10^ Mq). To some extent, the lack of higher mass BHs among NLSls in 
general may just be a 'selection' effect which arises from the standard definition of NLSls 
as objects with FWHMh/? < 2000 km/s. Since there is a clear correlation between H/3 
luminosity and H/3 width (e.g.. Fig. 15 of Veron-Cetty et al. 2001), higher luminosity AGN 
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should have broader Balmer hnes finally exceeding FWHM > 2000 km/s and thus would no 
longer be called NLSls. Ideally, any NLSl defining criteria should incorporate that effect 
(Wills et al. 2000). 

In any case, while we cannot exclude that NLSls with higher BH masses do exist, the 
fact that the radio-loud objects of our sample cluster in a previously almost unpopulated 

parameter space of Mbh and R is still remarkable. We caution there are indications that 
2-3 galaxies of the sample are actually beamed sources (Sect. 5.7.1). Any contribution 
of beamed emission to the optical band makes our estimate of black hole mass an upper 
limit, while, on the other hand, we underestimate the black hole mass if the BLR clouds 
are arranged in a disk rather than a sphere (Vestergaard et al. 2000). Among the sources 
of our sample, SDSSJ172206. 03-1-565451. 6 is the radio loudest with 'reliable' radio index 
measurement, since no indications for beaming and likely no absorption in the optical band 
(Komossa et al. 2006). 

5.4. Accretion rate 

With knowledge of BH mass and X-ray luminosity it is possible to perform an estimate 
of the accretion luminosity relative to the Eddington value. We intentionally use the X-ray 
luminosity here, and we do not apply any bolometric correction, since it is unlikely constant 
for all NLSl galaxies, and since, in particular, the EUV part of their SED is unknown^. 

We find that the ratio of (0.1-2.4 keV) luminosity to Eddington luminosity (Tab. 2) 
ranges between LjL^dd = 0.02 (RXJ23149+2243) and 0.7 (RXJ16290+4007) without any 
bolometric correction which might typically be an additional factor 5-10 (Tab. 14 and 15 of 
Elvis et al. 1994). An exception is SBS1517+520 with L^/LEdd = 0.007. We suspect that 
this very faint X-ray source is absorbed in X-rays by a (dust-free) absorber which would also 
explain its very fiat powerlaw with index ~ —1.1 which is formally obtained when fixing 
the absorption at the Galactic value. We excluded TEXll 11 1-^329 and 1RAS09426-M929 
from this estimate since no meaningful X-ray upper limit could be obtained. 

For comparison with Lboi/-^Edd values reported in the literature we also estimated bolo- 



^There have been suggestions that the X-rays connect in one 'big blue bump' to the optical band on 
the one hand, but findings that NLSls are under luminous in the UV on the other hand, and inspections of 
individual objects lead to expectations of flat EUV-SEDs in some cases, soft-excesses in other cases. Ideally, 
bolometric corrections would therefore be done for each object individually, based on a measurement of the 
multi-wavelength SED on the one hand, and detailed emission-line modeling on the other hand in order to 
constrain the EUV-part of the SED from emission-line ratios and H/3 line flux. 
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metric luminosities based on the optical luminosity at 5100A, and a constant bolometric 
correction factor of 9 (e.g., Sect. 2 of Warner et al. 2004, Greene et al. 2006); i.e. Lboi = 
9 X ALa(5100A). This gives an average Lboi/-Z^Edd=l-2 (again excluding TEXlllll+329 and 
IRAS09426+1929, plus 1RAS20181-2244) which is higher than the NLSl average of Warner 
et al. (2004; 0.67, but note that they used C1V1549 rather than B.(3 for BH mass estimates) 
but not as high as observed in the most luminous quasars (e.g., Warner et al. 2004, Shemmer 
et al. 2004). 

5.5. Radio-Fell connection, compcirison with Population A/B sources, and the 
radio- loud radio-quiet bimodality in 'Eigenvector 1' 

Among broad-line quasars the nature of the Fell emission has been much discussed in 
the literature with respect to radio properties (e.g., Joly 1991, Miller et al. 1993). Two 
models were considered to explain the connection between radio and Fell emission in radio- 
loud sources: an orientation-dependent model of Fell emission (e.g. Jackson & Browne 
1991), and a jet model in which Fell is produced when the jet interacts with the ambient 
medium (Joly 1991). Boroson & Green (1992) and Boroson (2002) find that radio loudness 
is statistically linked to the physical parameters that drive "Eigenvector 1" in the sense 
that quasars with strong [0111] and weak Fell are preferentially radio-loud while strong Fell 
emitters are generally radio-quiet objects (for an example of an exception, see Zhou et al. 
2002). 

The radio- loud NLSls of our sample do not fit in this scheme in the sense that their Fell 
emission does not show a clear dependence on radio loudness, but almost covers the whole 
observed range in NLSl galaxies (Fig. 3; cf. Fig. 4 of Marziani et al. 2001). In particular, 
RXJO 134-4258, the ~4th radio-loudest of our sample, is one of the strongest Fell emitters 
known. 

Sulentic et al. (2003) demonstrate that the radio-loud radio-quiet distinction of AGN is 
obvious in their 'Eigenvector 1' which links the width of the broad component of YLf] with the 
equivalent- width ratio of Fell to H/3 (their Fig. 3), in the sense that radio-loud and radio- 
quiet objects occupy different domain spaces. NLSl galaxies are not among their sample. 
Sulentic et al. define two populations of AGN according to their emission-line properties. 
'Population A' are sources with FWHMh/Ji, < 4000 km/s, while FWHMn/jb > 4000 km/s 
in 'population B' sources. In their sample, radio loudness is almost exclusively confined to 
population B. The radio-loud NLSl galaxies of our sample fall into the rarely populated 
regime of Fig. 3 of Sulentic et al., most of them in the 'gap' between radio- louds and 
radio-quiets, but not to an extent that the bimodality is obliterated. 
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5.6. Inclination 

Orientation plays an important role in AGN unification scenarios of both Seyfert galaxies 
(Antonucci 1993, Elvis 2000) and radio galaxies (Urry & Padovani 1995). Inclination also 
enters black hole mass and accretion rate estimates and consequently affects black hole mass 
- bulge mass relations of NLSls and their cosmological implications (see e.g. the discussion 
in CoUin & Kawaguchi 2004). 

The orientation of NLSls w.r.t our line of sight has repeatedly been addressed, and ar- 
guments and evidence in favor of and against a pole-on view have been presented. Originally, 
for instance, a pole-on view was suggested to explain the small width of the Balmer lines 
(e.g., Osterbrock and Pogge 1985, Puchnarewicz et al. 1992, Bian & Zhao 2004). The disk 
emission of Fell is also anisotropic and Fell emission is strongest when the disk is viewed 
facc-on (e.g., Miller et al. 1993, Marziani et al. 01). A pole-on orientation was also suggested 
for individual NLSls with extreme variability (e.g., PKS0558-504: Remillard et al. 1991). 
Correcting for these orientation effects, systematic or random in nature, is also important 
in the context of black hole mass estimates and black hole mass - bulge mass relations of 
NLSls, which make use of Balmer line widths. On the other hand. Smith et al. (2004) 
argued that, based on their polarization properties, there are no indications that NLSls are 
preferentially viewed face-on. 

Radio observations of NLSl galaxies allow us to re-address orientation issues. Among 
our sample, there are 2-3 NLSls which are hkely beamed (SDSSJ094857.3+002225, RXJ16290+4007, 
and possibly PKS0558-504; Sect. 5.8.1), implying a pole-on view onto the 'central engine'. 
If the width of the Balmer lines was dominated by orientation, we would then expect these 
objects to show very narrow Balmer lines. This is not the case (FWHMH/3b = 1710 and 
1800 km/s, respectively for SDSSJ094857.3+002225 and RXJ16290+4007), though. Larger 
samples of beamed NLSl galaxies are needed to put this on a firmer statistical basis. 

5.7. Frequency of radio-loud NLSl galaxies 

Of the 128 NLSl galaxies in the Veron quasar catalogue, 90% are located within the 
NVSS survey area. Among these, 7% are formally radio- loud (it!i.4^19), while only 2.5% 
exceed it! = 100. One or a few of these sources are most likely significantly absorbed and one 
or two are at the border between Seyfert 1 and NLSl galaxy. Given a fraction of ~13-20% 
of radio- louds among optically selected quasars (e.g., Visnovsky et al. 1992, LaFranca et al. 
1994, Falcke et al. 1996), radio- loud NLSls are more rare than radio-loud quasars. 

In order to compare directly the fraction of radio-louds among narrow-line and broad 
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line type I quasars listed in the VQC, we have also extracted the broad line quasars from 
the VQC, cross-correlated them with the NVSS catalogue, and determined the fraction of 
radio-louds the same way we did for the NLSl galaxies. We find that ~90% of the objects 
are located within the NVSS survey area. Among these, ~20% are radio-loud (i?i.4>19), and 
~14% exceed R = 100. Again, this indicates a lack of radio-louds among NLSls, particularly 
at high values of radio index R. 

When comparing the frequency of radio-loud NLSls with radio- loud type 1 quasars, 
two things have to be kept in mind. Firstly, the mere definition of NLSls as objects with 
FWHMh/3 < 2000 km/s likely biases the observed black hole mass distribution towards 
excluding the most massive black holes. Since there are indications that radio-loud AGN 
have, on average, higher BH masses (e.g., Metcalf & Maghochetti et al. 2006), we may 
automatically miss some of the most massive radio-loud NLSls because they escape the 
definition as NLSls and are rather called Seyfert Is. 

Secondly, specific to the present sample, great attention was paid to each single object 
in that cautious comments were issued in case of evidence for excess extinction, and in that 
line profiles were fit with a number of different models. The same detailed object-by-object 
analysis should be applied to comparison samples of quasars, when citing the fraction of 
radio-loud quasars in comparison with that of radio-loud NLSl galaxies. 

Keeping this in mind, it still appears that radio-loud NLSls are less abundant than 
radio-loud quasars, and we proceed in attempting to explain this distinction. 

The cause of radio loudness in quasars is still unknown, even though a number of different 
models have been suggested and investigated including the infiuence of spin (e.g., Wilson & 
Colbert 1995, Blandford 2000, Ye & Wang 2005), black hole mass (e.g., Laor 2000, Metcalf 
& Maghocchetti 2006), the rate of coohng gas in the host galaxy (Best et al. 2005) and 
beaming in radio-intermediate FSRQs which were suggested to be relativistically boosted 
radio-quiets (Miller et al. 1993, Falcke et al. 1996). 

Two key questions regarding the radio properties of the NLSl galaxies are: Why is the 
phenomenon of radio loudness rarer in NLSl galaxies than in quasars? And, related to that: 

Is the mechanism which drives the radio loudness the same in both types of objects, or is a 
fundamentally different mechanism at work in the few radio-loud NLSls? We discuss several 
models in turn. 
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5.7.1. Beaming 

We first explore the possibility that NLSls might be radio-quiet but intrinsically beamed 
sources. In this scenario, one then expects that other optical and X-ray properties of the 
radio- loud NLSls are affected as well. For instance, Fell emission might be extreme if it 
originates in a disk which is viewed face on (e.g.. Miller et al. 1993). Also, X-ray soft 
excess emission is expected to be stronger if the sources are viewed face on for the case of a 
geometrically thick accretion disk (Madau 1998). At the same time there might be a flatter 
non-thermal contribution to X-rays from the jet emission. 

While a few sources of our sample are indeed good candidates for beaming, others are 
not. RXJ16290+4007 and SDSSJ094857.3+002225 are likely beamed, but only RXJ16290+4007 
is radio- 'intermediate'. RXJ16290-I-4007 hkely is a blazar of FSRQ type (e.g., Padovani et 
al. 2002) and SDSSJ094857.3+002225 exhibits strong variability in the radio band (Zhou 
et al. 2003). These are the only two sources among our sample with very fiat radio spectra 
(Tab. 1). Apart from that, PKS0558-504 exhibits Timisual X-ray flaring activity and there 
have been repeated speculations that this may be related to beaming effects in a jet (Remil- 
lard et al. 1991, Gliozzi et al. 2001, Wang et al. 2001). An origin in the accretion disk 
is an alternative, though, and this latter model is preferred by Wang et al. (2001) and by 
Brinkmann et al. (2004) to explain the XMM data. 

Other sources of our sample are steep spectrum radio sources where strong beaming 
is generally not expected. X-ray spectra are steep (even those of the two FSRQs), arguing 
against the dominance of a non-thermal Doppler-enhanced contribution to the soft X-ray 
spectrum. Beamed sources more typically show Fx — —1.5.. — 1.6. 

5.7.2. Accretion mode 

While estimated Eddington accretion rates are in line with other findings that NLSl 
galaxies are accreting close to or above their Eddington hmit, the radio loudness of such 
sources is unexpected at first glance, given observations that radio emission is suppressed 
in accretion high-states in some types of Galactic and extragalactic sources. According to 
Maccarone et al. (2003, and references therein). Galactic X-ray binaries in soft/high-state 
and AGN show quenched radio emission for relatively high accretion rates. Even though the 
analogy between Galactic black holes and NLSl galaxies has been drawn repeatedly (e.g., 
Pounds et al. 1995), it is still unclear how far it reaches. Greene et al. (2006, and references 
therein) more generally suggested that radio loudness is anti-correlated with Lboi/-^Edd (their 
Fig. 3) in a sample of Seyferts and PG quasars, albeit with quite some scatter. Whether 
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this trend still holds for high-redshift, highly accreting quasars has yet to be explored. So 
far, these observations suggest, whatever the mechanism to suppress radio emission for high 
Lboi/LEdd could also be responsible for the lower fraction of radio-louds among NLSls. 

While it is generally agreed upon that jet formation (and thus radio emission) is inti- 
mately linked to the presence and structure of an accretion disk, the actual mechanism of 
launching a jet from the inner disk regions is not yet really understood (see, e.g., Celotti & 
Blandford 2001 and Meier 2003 for recent reviews). One may speculate that the observed 
differences originate from differences in the accretion mode between stellar mass and super- 
massive black hole accrctors and between different classes: It is generally agrcccd upon that 
at relatively low accretion rates (in terms of the Eddington rate), radiatively inefficient ac- 
cretion flows prevail (e.g., Beckert & Duschl 2002). Independent of yet not fully understood 
details of these flows, they seem to be potentially efficient in producing radio emission them- 
selves (e.g., Yi & Boughn 1998, Narayan et al. 1998). At rates closer to, but still well below 
the Eddington limit disks turn optically thick and geometrically thin (sec Wang et al. 2002 
for a suggestion how radio-emission might be enhanced in this case). Barring details at these 
low accretion rates the properties of the flows seem to be rather similar between stellar mass 
and super-massive black hole accretors. 

When the accretion rate reaches or even surpasses the Eddington rate, however, the 
physical state of the accretion process becomes less clear. One has to keep in mind that 
the Eddington rate has been derived for a perfectly spherically symmetric situation. The 
situation in accretion disks, however is rather different from such an ideal case. While at 
rates close to the (spherical) Eddington limit, accretion disks in all likelihood are no longer 
geometrically thin {slim disks, see, e.g., Collin & Kawaguchi 2004), their deviation from 
sphericity is even more important. This can lead to simple geometrical effects allowing for 
super- Eddington accretion rates (Heinzeller and Duschl, in prep.), but also to instabilities 
in the structure of the disks, for instance the photon bubble instability as discussed by 
Begelman (2002), and in all likelihood to combinations of such geometrical and physical 
effects (e.g., Okuda et al. 2005). 

An additional potentially relevant difference between stellar mass and super-massive 
black hole accretors may be the mass of the disk. While in the stellar case the masses of the 
accretion disk are certainly well below that of the accreting body, this is considerably less 
clear for the AGN case, and may be even different for different AGN types. In the innermost 
disk regions, selfgravity certainly plays no role. Already at moderately large radii, however, 
this may change. With such changes, long-term modifications of the mass supply rate for 
the inner regions comes along. Different levels of selfgravity in the disks may very well lead 
to differences in the accretion process. 
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Therefore, different accretion modes have the potential to explain the observed differ- 
ences in radio loudness between Galactic and AGN sources on the one side, and among 
the different AGN classes on the other side. Before one, however, can unambiguously iden- 
tify whether, and if so, which aspect of the accretion mode gives rise to the differences, a 
much better understanding is required of how these different accretion modes relate to jet 
formation. 



5. 7. 3. Black hole spin 

Another parameter which could account for radio loudness but at the same time leave 
optical emission-line parameters and correlations mostly unaffected, is black hole spin. It is 
expected to have a strong direct influence on radio jet emission if the radio jet is powered by 
the extraction of black hole rotational energy (Blandford & Znajek 1977), or may indirectly 
be relevant if the jet is powered by the rotational energy of the inner disk (Blandford & 
Payne 1982). 

If, indeed, rapidly spinning black holes arc responsible for radio loudness of NLSls and 
quasars, why then would NLSls on average harbor less rapidly spinning black holes than 
Seyferts (to account for the lower frequency of radio-loudness among NLSls)? Both BH-BH 
mergers and accretion will affect BH spin (e.g., Merritt 2002, Hughes & Blandford 2003, 
Volonteri et al. 2005). It has been suggested that NLSl galaxies are less evolved than 
Seyferts in that their black holes are still in the process of growing by rapid accretion (e.g., 
Mathur 2001, Grupe & Mathur 2004, Botte et al. 2004, Mathur & Grupe 2005). Since 
accretion processes are known to be very efficient in spinning up the accreting black holes 
(Volonteri et al. 2005, and references therein) one then indeed expects that the flnal stages, 
'normal' Seyferts and quasars, should have, on average, more rapidly spinning black holes 
than the NLSls. 

At the same time, spin cannot be the only quantity making objects radio-loud, else all 
quasars must be more radio-loud than NLSls. Thus, spin would be a necessary but not 
sufficient condition. If, indeed, spin plays a role, then the few radio- loud NLSls should be 
those with the more rapidly spinning BHs, already more evolved than others. One may thus 
expect that they he closer to the black-hole mass - velocity dispersion (Mbh — c) relation of 
Seyfert galaxies than other NLSls. We have plotted the objects of our sample in the Mbh — c 
plane of Mathur & Grupe (2004; their Fig. 1) and indeed find that they are systematically 
closer to, or on, the Mbh — c relation followed by Seyferts. 

Future X-ray observations of the iron line may help to search for signs of more rapidly 
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spinning BHs in the radio-loud NLSls, but less rapidly spinning BHs in the NLSl population 
as a whole. 

In summary, while we cannot fully exclude beaming to operate in all sources, the source 
properties are generally such that this mechanism is not favored (with the three exceptions 
mentioned above). Accretion mode, possibly in combination with spin, is a possibility. 

5.8. Future work 

The results of this study lead us to identify a number of important follow-up observa- 
tions: 

High quality optical spectra of the 11 radio-loudest objects of our sample do not yet 
exist, or are not yet published. Spectra of high resolution and S/N would allow a search for 
similarities/differences in emission lines between radio-loud and radio-quiet NLSls the same 
way it is now rigorously done for broad- and narrow-line objects (e.g., Veron et al. 2001, 
Sulentic et al. 2002, Bachev et al. 2004). Especially important is a rehable decomposition 
of the faint Balmer lines like H7 and US into broad and narrow component, and dcblending 
from [OIII]4363, such that extinction measurements become possible for those objects which 
have Ha redshiftcd out of the observable optical band. This knowledge is essential for the 
assessment of radio loudness (in reddened objects, radio loudness is overpredicted) , and for 
black hole mass calculation. 

Good optical spectra of the radio-loud and almost radio-loud NLSls, and a comparison 
sample of quasars, will enable BH mass and bulge mass estimates using H/5 and [OIII] or 
low-ionization lines, and ideally also absorption features (e.g., Grupe & Mathur 2004, Botte 
et al. 2004, Greene & Ho 2005) which will then allow us to test more rigorously whether 
claimed systematic trends of higher black hole masses in radio-loud vs. radio-quiet quasars 
extend to the NLSl population. 

X-ray observations of all sources are needed in order to measure accurate spectral shapes, 
few of them available so far, and to determine the amount of absorption. Deeper X-ray 
observations are required to measure iron-line profiles in search for (less) rapidly spinning 
black holes in NLSl galaxies. 

How much of the scatter in radio spectral index of the NLSls is real and how much 
is caused by the non-simultaneity of the data, should be tested with simultaneous radio 
observations at different frequencies and monitoring for variability. The latter will also 
constrain the number of beamed sources. Radio observations of higher spatial resolution 
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should be carried out in order to confirm the compactness of the NLSls of our sample and 
thus also to check how far the similarity with known CSS sources goes. 

Larger samples of AGN with a 'relaxed' NLSl classification criterion where the cut 
in FWHM is dependent on luminosity will enable us to address the question whether the 
fraction of radio- loud NLSls systematically changes (increases) with FWHMh^ and whether 
(radio-loud) NLSls with higher BH masses emerge. The largest sample to date of NLSls 
from the SDSS data base, but still using an FWHM cut-off, is presently being compiled by 
Zhou et al. (priv. comm.). 

Apart from future multi- wavelength follow-up observations of the small sample of radio- 
loud NLSls presented here, Virtual Observatory tools will become more and more important 
to study trends throughout the whole NLSl and AGN population, since they will allow the 
extraction of multi-wavelength information from existing data bases in a homogeneous way 
which then enables us to study bulk properties of populations on the one hand, and to 
identify exceptional individual sources on the other hand. 

Comparison with the work of Whalcn et al. (2006): While our paper was being refereed, 
a paper by Whalen et al. (2006) appeared on the astro-ph preprint server. These authors 
present optical follow-up spectroscopy of FIRST radio sources identified as NLSl galaxies, 
including 16 candidate radio-loud NLSls (>50% of their sample has radio powers < 10^^ 
W/Hz). Their findings are consistent with ours in that both works do not find strong 
correlations between optical line parameters and radio-loudness, and black holes masses 
show similar mass ranges. Being a radio-selected sample, the fraction of radio-loud NLSl 
galaxies is higher in the sample of Whalen et al.; on the other hand a number of the sources of 
Whalen et al. are below a radio power of 10^^ W/Hz, and we note that sources significantly 
above a radio index i?=1000 are missing in both studies. 



6. Summciry and conclusions 

We have presented a study of the radio properties of NLSl galaxies, with emphasis on 
the search for radio-loud NLSls. The results of this study can be summarized as follows: 

11 candidate radio- loud NLSl galaxies were identified, significantly increasing the num- 
ber of known radio-loud NLSls. These are radio loud judged on both, radio index and radio 
power. Radio- detected NLSls cover several orders of magnitude in radio index R. The 
distribution is rather smooth below R ~ 10. 

The fraction of detected radio- loud NLSls is ~ 7%, smaller than the fraction of radio- 
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louds among quasars. Only 2-3% of the NLSls exceed a radio index R > 100 (SDSSJ172206.03+565451.6, 
SDSSJ094857.3+002225 and possibly the absorbed source TEXlllll+329). This could 
partly, but likely not fully, be a 'selection effect' arising from the NLSl definition as objects 
with FWHMh/3 < 2000 km/s independent of luminosity. 

The radio-loud NLSls are generally steep spectrum sources and compact. Two have 
inverted spectra. NVSS and FIRST radio fluxes are consistent with each other in most 
cases, indicating little variability. 

Accretion rate estimates show that most sources likely accrete close to or above the 
Eddington limit. Black hole masses are generally at the upper end observed in NLSl galaxies 
(Mbh ~ (2 — 10) 10'' M0, for objects with optical spectra of good quality), but still small 
when compared to some samples of radio-loud quasars. In particular, the NLSls of our 
sample fill up a previously scarcely populated area of Mbh — R diagrams. 

Optical properties of the radio-loud NLSls are similar to the NLSl population as a 
whole. In particular, their Fell and [OIII] emission covers almost the whole range observed 
in NLSl galaxies. Their radio properties, however, extend the range of radio- loud objects to 
those with small FWHMh/j. The fraction of [OIII] 'blue-winglers' is high. X-ray properties 
are consistent with the NLSl population as a whole, with Fgx — —2.. —3.5 and Lgx — 10^^~^^ 
erg/s; some sources highly variable, others constant throughout the observation interval. 

The starburst contribution to the radio emission, estimated for those sources with IRAS 
60/U detections, is small or completely negligible. The AGN is expected to dominate any 
starburst in radio emission by factors 5-115. 

We do not find indications for beaming in most of the sources, arguing against the 
explanation that we see intrinsically radio quiet but beamed sources. Accretion mode and 
spin might be factors to account for the lower frequency of radio-loud NLSl galaxies than 
quasars. 

Future studies of radio- loud NLSls and follow-up observations of the sample presented 
here will provide important information on distinguishing between different NLSl models, 
on (the existence of) the radio-loud radio-quiet dichotomy of AGN, and on accretion disk - 
jet models. 
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Fig. 1. — Radio-detected NLSl galaxies selected from the VQC, ordered according to radio loud- 
ness. The solid line marks the "deviding line" between radio-loud and radio-quiet objects, Ria = 19 
(corresponding to i? = 10; see text for definitions). Calculation of Ria is generally based on USNO- 
Bl blue magnitudes mB2 (squares) and FIRST (plotted in black) and/or NVSS (plotted in blue) 
radio detections, except RXJ 0134-4258 (PMN data at 4.85 GHz, plotted in red), PKS0558-504 
(PMN data, red; and SUMSS data at 0.8 GHz, green), and RXJ22179-5941 (SUMSS data, green). 
Several FIRST/NVSS objects have radio detections at other wavelengths. These are not over- 
plotted here, but are used for estimates of radio spectral indices (Tab. 1). In order to have an 
estimate on the uncertainty in me - either due to real source variability or measurement uncer- 
tainties - we also show i?i.4 calculated using mei (circles) of the USNO-Bl catalogue, me of the 
USN0-A2 catalogue (triangles), mBj taken from the GSC2.2 catalogue (crosses), and corrected 
4400A fluxes obtained from the SDSS, ESO and Xinglong spectra (small filled circles). The dotted 
line indicates a change in blue magnitude by Ame = ±1 mag. Plotted ratios are for the following galaxies, 
from left to right: N=l: TEXlllll+329, 2 = SDSSJ094857.3+002225, 3 = SDSSJ172206. 03+565451. 6, 4 = IRAS20181-2244, 
5 = RXJ0134-4258, 6 = PKS0558-504, 7 = IRAS09426+1929, 8 = SBS1517+520, 9 = RXJ16290+4007, 10 = 2E1346+2637, 
11 = RXJ23149+2243, 12 = IRAS11598-0112, 13 = RXJ22179-5941, 14 = IRAS00275-2859, 15 = IRAS20520-2329, 16 = 
RXJ01354-0426, 17 = MS12510-0031, 18 = IRAS11058+7159, 19 = HS0710+3825, 20 = SBS1152+523, 21 = SBS1126+516, 
22 = IRAS13349+2438, 23 = RXJ15308+2026, 24 = RXJ16196+2543, 25 = RXJ11425+2503, 26 = RXJ15475+1024, 27 = 
FBS1002+437, 28 = PG1404+226, 29 = RXJ17025+3247, 30 = RXJ12257+2055, 31 = PG1448+273, 32 = MRK478, 33 = 
PHL1811, 34 = IZwl, 35 = IRAS13224-3809. 
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Fig. 2. — Distribution of the radio index i?i 4 of the NLSl galaxies of our sample. i?i 4 was 
calculated using NVSS radio fluxes and USNO-Bl mB2 magnitudes, except if only FIRST 
fluxes and/or mei magnitudes were available which were then used. The horizontal fat 
bar marks the range within the radio-loud regime which is sometimes referred to as 'radio 
intermediate'. 
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Fig. 3. — Distribution of the radio-loud NLSl galaxies of our sample (filled squares) in the 

FeII/H/9toti-[0III]/H/?toti diagram of NLSl galaxies (open squares; Xu et al. 2006). Square 
size codes radio-loudness; the larger the square diameter the radio-louder the source. The 
source IRAS20181-2243 is off the plot, indicated by the arrow. 
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Fig. 4. — The dependence of radio-loudness on black holes masses for the radio-loud NLSl 
galaxies of our sample for which spectra were available to us (fat circles) and based on spectra 
in the hterature (open circles) . The dashed and sohd hues mark areas populated by the bulk 
of the radio sources of the samples of Laor (2000; solid line) and Lacy et al. (2001, dashed 
line). The open double-circle marks TEXlllll-l-329 which is possibly strongly absorbed. 
See text for further details, including comments on individual sources. If both FRIST and 
NVSS radio observations were available, radio indices shown in this plot were calculated 
based on the NVSS data. 
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Table 1: Summary of the radio properties of radio-loud NLSls and candidates, plus PKS 2004- 
447 (a peculiar NLSl, Oshlack et al. 2001). The range in radio index R of each galaxy is 
defined by the choice of the blue magnitudes from different epochs/catalogues as in Fig. 1, 
and is based on FIRST observations if available, z is the redshift, z/r gives the observed 
frequency at which R was calculated, is the radio powerlaw spectral index {fi, oc z^"^"), 
and 1^1,1^2 specify the frequency interval in which was calculated. 



name 




z 


R 


GHz 






1^1,1^2 

GHz 


comments, references 


TEXlllll+329 


0. 


,189 


1249. 


,.445 


1.4 


-0. 


,56 


0.33,1.4 


optically absorbed 














-1, 


,24 


1.4,4.85 
















< -2. 


,07 


1.4,15 


upper limit at 15 GHz ([1]) 


SDSSJ094857+002225 


0. 


,584 


793., 


,194 


1.4 
















>1000 


5 


0, 


,59 


2.7,4.85 


[2], variable 


SDSSJ172206+565451 


0. 


,425 


773. 


,.70 


1.4 


-0. 


,69 


0.33,1.4 


[4] 


RXJ0134-4258 


0. 


,237 


178. 


.36 


4.85 


-1. 


,43 


4.85,8.4 


see also [3] 


IRAS20181-2244 


0. 


,185 


52., 


.22 


1.4 


-0. 


,51 


0.35,1.4 




RXJ16290+4007 


0. 


,272 


41., 


.21 


1.4 


0, 


.42 


1.4,4.85 










68., 


.35 


4.85 










IRAS09426+1929 


0, 


,284 


28., 


.16 


1.4 










PKS0558-504 


0. 


,137 


16., 


.14 


4.85 


-0, 


,45 


2.7,4.85 
















-0. 


,38 


0.84,2.7 




RXJ23149+2243 


0. 


,168 


18. 


,.8 


1.4 










2E1346+2637 


0. 


,915 


18. 


,.6 


1.4 










SBS1517+520 


0. 


,371 


27., 


.12 


1.4 










PKS 2004-447 


0.24 


6320., 


.1710 


4.85 


-0.67 


ATCA 


[5], variable 



References: [1]: Nagar et al. (2003), [2]: Zhou et al. (2003), [3]: Grupe et al. 2000, [4]: 
Komossa et al. 2005, [5]: Oshlack et al. 2001 
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Table 2: Summary of optical and X-ray properties of the radio-loud NLSls of our sample, wh/? 
and W[oiii] correspond to FWHMh/3 and FWHMjomjsoo?) respectively, in km/s. The three 
entries for FWHMh/3 arc based on a single-component Gauss fit, a direct fit, and a two- 
component fit with narrow component fixed to the same width as [Olll]5007, respectively. 
[Olll]5007, Fell4570 and Ha are intensity ratios relative to H/?g. 'CR' and Fx refer to the 
ROSAT PSPC countrate in cts/s and the soft X-ray photon index, respectively. Column 
'ref provides the references for the optical spectroscopy. Entries of the first six objects are 
based on our (re-) analysis of SDSS and other spectra. The other entries were taken from 
the literature. In the latter case, multi-component fits to the Balmer lines are generally not 
available, and only results from single Gaussian fits to the emission lines are listed. 



name 


Wh/3 


W[oiII] 


[OIII] 


Fell 


Ha 


CR 


Fx 


ref 


SDSSJ094857+002225 


1420/ 960/1710+wJqjj] 


400^ 


0.2 


1.3 




0.04 


-2.2 


[1] 


SDSSJ172206+565451 


1580/1490/1990+W[oiii] 


490 


0.6 


0.7 




0.08 


-2.. - 3 


[1],[6] 


RXJ0134-4258 


930/1040/1860+W[om] 


510 


0.04 


3.2 


3.4 


0.20 


-2.2^ 


[1] 


RXJ16290+4007 


1260/1290/1800-Fw[oiii] 


310 


0.4 


1.1 


3.3 


0.83 


-3.0 


[1] 


RXJ23149-F2243 


1670/1630/2680-Fw[om] 


620^ 


0.2 


1.4 




0.11 


-1.9 


[1] 


SBS1517+520 


3220 /2030 /3810-Fwfoiiii 


1140 


1.1 


0.5 




0.015 


-1.1 


[1] 


TEXlllll+329 


1980 


1480 


0.75 


1.1 




<0.04 




[2] 


IRAS09426+1929 


<2000 










< 0.03 






PKS0558-504 


1500 




0.1 


1.56* 




4.20 


-2.8 


[3] 


2E1346-F2637 


1840 


1130^ 




1.0 




0.08 


-3.4 


[4] 


IRAS 20181-2244 


460 


695 


3.4 


0.7 


6.7 


0.04 


-2.6 


[5] 



^ Width of narrow component fixed to [Oil] 3227 since [OIII] 5007 is too weak; ^the reported 
FWHM is that of [Oil] rather than [OIII] (in case of 2E1346+2637 the [Oil] width is not 
de-convolved from the instrumental profile). ^The [OIII]profile has a very strong and broad 
blue wing. The FWHM reported here is that of the narrow core of the line only. ^During 
the RASS, the spectrum was significantly steeper (Fx = —4.3). 



*Sum of both Fell complexes. 

References: [1]: this paper, [2]: Zheng et al. 2002, [3]: Remillard et al. 1986, [4]: 
Puchnarewicz et al. 1996, [5]: Halpern & Moran 1998, [6]: Komossa et al. 2005 
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Table 3: Estimated black hole masses and Eddington ratios L^/L^^d (multiplied by a factor 
10) for NLSls with SDSS or other spectra (broad component of H/3 used for BH mass 
estimate; first 6 objects) and other NLSls (single Gauss fits to H/5 from literature used for 

BH mass estiinates). 
name -Mbh lOxLx/LEdd comments 



Mo 


SDSSJ094857+002225 


410^ 


4.5 


source might be beamed; see also [1] 


SDSSJ172206+565451 


310"^ 


3 




RXJ0134-4258 


410^ 


1 




RXJ16290+4007 


210^ 


7 


source might be beamed 


SBS1517+520 


110^ 


0.07 


Ukely X-ray absorbed Lx/L^dd is lower limit 


RXJ2314+2243 


810^ 


0.2 




PKS0558-504 


910^ 


7 


see also [2] 


TEXlllll+329 


210"^ 




likely optically absorbed 


IRAS09426+1929 


< 910^ 




FWHMh/3<2000 km/s assumed; probably optically absorbed 


2E1346+2637 


810^ 


4 




IRAS 20181-2244 


310^ 


6 





[1]: Zhou et al. (2003), [2]: Wang et al. 2001 



